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INVARIANT  IMBEDDING  AND  GENERALIZED  TRANSPORT  THEORY- - 
A  BASIC  STOCHASTIC  FUNCTIONAL  EQUATION 

Richard  Bsllasui  and  Robert  Kalaba 
The  RAND  Corporation.  Santa  Monica ,  California 

and 

G.  Milton  Wing,  Lo*  Alamo*  Scientific  Laboratory 

Loa  Alamo* ,  Nev  Mexico* 

1.  Introduction .  In  a  eerie*  of  recent  paper*,  the  principle  of 

Invariant  Imbedding  ha*  been  utilised  in  the  *tudy  of  a  variety  of 

12  3  l(  f 

physical  prooesee*:  radiative  tran*fer,  ’  (*ee  Prei*endorfer  ’ "j 

neutron  tianeport,^  ^  random  walk  and  scattering .  and  vae 
propagation.11 

The  ala  of  the  preeent  paper  1*  to  extend  previous  result*  and 
technique*  *o  a*  to  Include  an  extensive  category  of  transport 
processes  involving  both  deterministic  and  stochastic  interaction 
general  gecmmtrles,  and  the  determination  of  characteristic  functions 
.nd  probabilities  as  veil  as  fluxes,  vhlch  Is  tc  say,  expected  values 

Utilising  the  principle  of  Invariant  Imbedding/  a  basic  stochastic 
functional  equation  vlll  be  derived.  vrcm  this  equation  appropriately 
specialised,  can  be  obtained  all  the  relations  pertaining  to  fluxes 
contained  in  the  foregoing  papers,  and  In  addition  corresponding 
relatione  for  characteristic  functions.  From  these  equations  for 
higher  moments  can  be  obtained.  As  far  as  ve  knew,  this  use  of 
stochastic  functional  equations  is  nev. 

For  intuitive  purposes,  the  reader  nay  use  as  a  model  neutron 
transport  theory.  In  order  to  emphasise  the  cosenon  features  of  a 
□umber  of  particle  processes,  ve  have  encased  the  problem  In  a  more 

abstract  setting.  In  this  vay  ve  obtain  an  extensive  generalization 

12 

of  the  fundamental  invariance  principles  of  Ambarxualan  and 
Chsuidrasekhar . 13 
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2.  Description  of  a  Generalised  Transfer  Process.  By  a  'particle'  we 
•hall  Man  a  state  vector  p  specified  by  a  position  coordinate  x.  The 
state  vector  contains  information  regarding  energy,  direction  of  motion, 
and  other  Information  required  to  specify  the  type  of  particle.  As 
the  particle  proceeds  through  a  me dint .  ip  engages  in  interaction  of 
de termini stia  or  stochastic  nature  vita  the  medium.  This  interaction 
leaves  the  medliai  unaffected,  but  is  equivalent  to  a  sequence  of 
transformations  of  the  state  vector,  and  in  sene  cases  of  the  number 
of  particles.  We  shall  assume  that  there  are  no  particle -particle 
interactions.  Although  me  ha  vs  begun  in  another  publication0  the 
study  of  collision  processes  using  invariance  principles ,  a  general 
formulation  of  these  more  complicated  processes  la  left  for  another 
time. 

In  order  to  use  functional  equation  techniques .  ve  utl lire  the 
concept  of  a  stratified  medli*.  In  a  finite -dimensional  space, 
conceive  of  a  family  of  surface# ,  of  one  leaa  dimension,  each  of 
which  is  specified  by  the  parameter  x.  TV. is  family  is  assumed  to 
have  the  property  that  its  members  can  be  used  to  partition  the  whole 
space  into  a  denimerable  set  of  strata  with  the  following  continuity 
property.  Each  atratv*  possesses  the  property  that  the  probability 
of  an  interaction  between  a  particle  and  the  medium  within  this 
strati*  can  be  made  arbitrarily  small  by  choosing  the  strati*  to  be 
bounded  by  arbitrarily  closely  neighboring  members  of  the  family  of 
surfaces.  Schematically . 


Pig.  1 
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We  shall  call  the  region  to  the  right  of  an  x- surface ,  an  x- region 

The  moat  common  stratifications  are  those  aecoapliahed  by  planes, 
spheres  and  cylinders,  but  more  general  families  of  surfaces  can 
equally  veil  be  used. 

Given  a  region  of  space  bounded  by  an  x-aurfaoe  and  an  0- surface , 
ve  wish  to  determine  the  maaber  and  nature  of  the  particles  existing 
to  the  left  of  the  x-surface.  'reflected'  particles,  and  those  existing 
to  the  right  of  the  0-surface,  ' transmitted '  particles,  as  a  result  of 
a  particle  specified  by  p  incident  on  the  left  as  in  Fig.  1. 

Let  us  now  introduce  some  random  variables  vhich  ve  shall  use  to 
describe  this  transport  prooess  analytically 


t(p,q;x)  -  the  random  number  of  particles  in  state  q 

reflected  from  the  x-region  over  all  .ime ,  due 

l  o  i) 

to  an  initial  particle  in  state  p  impinging 
upon  the  x-surface  from  the  left  at  time  rero 

Baa  direction,  velocity,  and  position  on  the  x-surface  are  all  contained 
vl thin  the  state  vector . 


r(p;x) 


1,  if  the  particle  in  state 
interaction  in  the  stratum 


stratim 


c-A,x]  . 


0,  otherwise  . 


p  is  involved  in  an 
x ,x-d  J  ,  or  in  the 


(2.2) 


•  (Pi -  1»  If  the  result  of  an  interaction  is  to 

produoe  k  particles  in  states 

*1  ’ **2  ’  ’  *  ‘ \  ’  *  "  1  ’  ^ . 

-  0,  otherwise 

By  the  stochastic  variables  *(1\p,q;x),  r^(p.z),  s  1/(p;q1.q?  ,.q^;x) 

ve  shal 1  mean  respectively  any  of  a  denumerable  set  of  variables  vith 
the  properties  described  above . 

Let  T(p,x)  denote  the  deterministic  change  in  state  caused  by 
passage  through  the  [z,x-d]  strati®. 

Finally,  ve  make  the  simplifying  assumption  that  there  are  only 
a  finite  set  of  possible  states.  TTils  permits  us  to  deal  vith 
characteristic  functions  rather  than  characteristic  functionals, 
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and,  in  any  cam,  la  the  type  of  aasiMptlon  required  to  carry  through 
any  computation* . 

3.  Verbal  Description  of  Process .  In  order  to  obtain  the  functional 
equation  of  the  following  auction  we  vlev  the  procee*  In  the  following 
faahlon.  A  particle  incident  upon  the  x- surface  undergoes  a  determ. n 
istlc  transformation  and  a  rtochartic  tranaf omation  in  the  Lx,x-A^ 
atratusi.  The  raault  of  the  date  mini  ati  c  tram*  format  1  on  la  to  con  ert 
it  from  a  state  p  to  a  atate  T(p,x).  The  result  of  the  stochaatlc 
tranaf ormat ion  la  to  produoe  a  random  number  of  p&rtlcies  In  a  random 
aet  of  atataa.  Each  of  theae  particles  acta  independently  of  the  others 
as  a  particle  incident  upon  the  (x-A)  surface  Aa  a  result  of  a 
particle  incident  upon  the  (x-A)  surface,  there  is  produced  a  random 
number  of  emergent  particles  In  a  random  aet  of  states  Each  of 
theae  undergoes  a  deterministic  and  stochaatlc  transformation  as  It 
travel*  through  the  j^x-A  xj  etratim  and  so  on 

In  deriving  the  baala  equation*,  however,  we  need  not  pay 
attention  to  particle*  which  are  the  result  of  more  thAn  one  determin¬ 
istic  or  stochaatlc  Interaction,  since  in  either  came  the  end  effects 
2 

are  of  order  A  .  It  Is  this  property  which  ocables  ua  to  derive 
relatively  simple  equation*  describing  quite  complicated  processes  , 
and  It  is  for  this  reason  that  we  employ  the  atratif icatlon  described 
above . 

Basic  Btochaatl c  Functional  Equation .  Using  the  invariant 
Imbedding  technique  and  enumerating  events  we  are  led  to  the  following 
stoehartlc  functional  relation: 
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1 


f  k  (1) 

•(p,qjx)  -  r(jjx)  ^•(p;q1,ci2..  •  •  ,\;*)  ,  Y2* 


-A) 


i-1 


i-n 


.(i) 


♦  (l-r(p;x))  £Z(l-rv  (q,x))  (n  -  *(T(p.x)fl.x-A) ) 
1-1 


1  — n 1 

♦  (l-r(p;x))  22  12 

M  '-1 


r  ■i^(qm;*)J  •^^(qm;w1  ,w2<  .  .v^jx) 


■  (Vi) 


(«•  1) 


,U)/„ 


.  22*  \v4-T"  (q)  .x-A) 

j  j 


♦  y(p.q;x) , 


(n*  -  t(T(p.x)^'x) ) 


vhere  y(p,q]x)  Is  a  contribution  f ra*  events  that  ha-e  probability 

2 

0(A  )  as  A — *0.  A  similar  equation  oan  be  derived  for  the  neutron* 

transmitted  through  the  stratus  [x,o]  . 


5.  Dlscu**lon.  Taking  expected  value*,  ve  derive  the  1  lux  equation* 
of  the  type  appearing  in  the  papers  cited  above.  Taking  the  expected 
value  of  (  ve  obtain  after  a  certain  amount  of  analytic 

manipulation,  0  or  re*  ponding  functional  equation*  for  the  characteristic 
function*.  Tbs  full  results  vill  be  presented  subsequently. 
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